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ABSTRACT 


The  effects  of  various  geometrical  boundary  perturbations 
on  finite-amplitude  acoustical  standing  waves  in  a rectangular, 
rigid-walled  cavity  were  investigated  using  non-linear  theory. 
The  standing  waves  that  exist  in  an  ideal  cavity  must  be 
corrected  when  the  boundaries  are  irregular .Three  specific 
examples  (stepped,  linear  and  wedged  perturbations)  were 
worked  out  to  demonstrate  the  corrections  (in  first  order) 
near  degeneracies  for  small  perturbations .Those  specific 
examples  were  compared  to  the  experiments . The  present 
theoretical  model  qualitatively  predicts  the  effect  of  the 
perturbations  on  the  behavior  of  the  nonlinearly  generated 
second  harmonic .However,  there  are  unexplained  quantitative 
discrepancies  between  experiment  and  theory  for  a couple  of 
cases . 
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1.  INTRODUCTION 


The  purpose  of  this  research  was  to  investigate  some 
of  the  effects  of  boundary  wall  perturbations  on  finite 
amplitude  standing  waves  in  a rigid-walled  rectangular 
cavity. The  investigation  was  prompted  by  an  examination 
of  the  experimental  results  of  Coppens  and  Sanders[33, 
the  research  of  DeVall[5]  and  of  Kilmer[43 , which  suggested 
the  existence  of  the  excitation  of  modes  other  than 
those  belonging  to  the  family  of  the  driven  mode. 
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2 . BACKGROUND 


A plane  elastic  wave  travelling  in  a non-dissipative 
fluid  will  change  waveform  as  predicted  by  the  relevant 
hydrodynamic  equations  [6],  [7]  .If  the  problem  is  extended 
to  absorptive  media, only  waves  of  relatively  high  amplitude 
will  change  waveforms  appreciably. 

At  the  Naval  Postgraduate  School ,Coppens  and  Sanders 
[3]  , Kilmer  [4]  ,and  DeVall  [5]  have  dealt  with  the  study 
of  finite  amplitude  waves  in  rigid-walled  rectangular  cavities. 

One  interesting  result  of  these  cavity  experiments  was 
the  appearence  of  excitations  of  modes  which  were  not 
family  members  of  the  driven  mode. For  example , assume  a.  rigid 
cavity  of  dimensions  Lx,Ly,Lz  is  driven  acoustically  at 
frequency  w,the  resultant  pressure  standing  wave  is  of  .the 
form 


coskxx  coskyy  coskzz  cos(wt+8)  (2.1) 

where  kx=NTT/Lx  , k =MTT/L  , kz  = LTT/Lz  (2.2) 

and  N,M,L  are  integers . Eq .( 2 . 1)  can  be  represented  by  the 
notational  shorthand 
(N,M,L/w,0) 

If  the  cavity  is  driven  to  excite  the  (0,M,0)  mode, then  the 
family  of  standing  waves  consist  of  all  of  those  of  the  form 
(0 ,nM,0/nw,9n) when  nw=nwQ  M Q. 

As  it  is  stated  in  £3]  ,"The  standing  waves  which  can  be 
excited  in  any  real  cavity  deviate  from  the  predictions  of 
the  linear  wave  equation  with  ideal  boundary  conditions 
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"for  the  following  reasons: 

(a) . The  presence  of  boundary-layer  losses  at  the  cavity 
surfaces  yields  a dispersive  contribution  to  the  wave  equation, 

(b)  . Geometrical  irregularities  alter  the  effective 
dimensions  of  the  cavity. 

Both  of  these  mechanism  can  be  treated  as  equivalent  as 
long  as  the  shift  in  frequency  are  so  small  that  the  actual 
resonances  are  close  to  the  theoretical  values  resulting 
from  the  classical  model."  These  are  treated  by  assuming 
the  dimensions  are  exact, and  the  apparent  phase  speed  is 
determined  on  that  basis. 

The  resonance  frequency  for  each  standing  wave  is 
defined  as  (0 

“n=c„I>x  kx>2  * <ny  V’  * <"z  (2-3) 

where  k's  are  given  by  Eq.(2.2)  and  n is  a shorthand  for 

the  set  (n.n „,n_)  where  n.n,  ,,n  are  integers , and  C is 
x y z x y z n 

the  apparent  phase  speed  appropriate  for  that  frequency. 

The  non-linear  wave  equation  applicable  to  this 
problem  can  be  obtained  as  follows: 

The  continuity  equation  for  wave  propagation  in  Eulerian 
coordinates  is 

*(/*)  + ! f=°  (2.4) 

this  equation  can  be  written  in  terms  of  the  condensation, 

”/o  )/ /o  as 

V.  C(  1+s)  uj  +|f  =0 


(2.5) 


The  equation  of  motion  in  Eulerian  coordinates  for  a 
contained  viscous  fluid  is 
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/>  [ 4 ( U . 7 ) a J * - * P * b V ( 7.  U ) - 7 V X v X U + 0 0, A T 


3t 

where 


P*  4£L{jLfbJt£[\  + *s+*-^r1  *] 


C2,ai 


(2.7) 


0DAT=0ther  dispersive  and  absorptive  terms  arising  from 
boundary  effects. 

Eq.(2.6)  can  be  rearrange  in  the  form  of 

|£+(ui.7)a  = 7-  ^ (2.8) 

where  the  operator  £ describes  those  physical  processes 
leading  to  absorption  and  dispersion. 

One  can  write  VxU=0  and  therefore  UL*  7$  , where  (J)  is  the 
velocity  potential , based  on  the  irrotational  velocity 
assumption .Hence , £ ll*  tv$  .Replacing  Us  7$  and  using  the 
condensation  ,Eq . (2 . 8 ).  can  be  written  as 

£v$+iv(vi)\j£vO+s)r''=  (2  9) 

Now, with  the  help  of  Eq.(2.5)  and (2. 9),  ard  after  a good 
deal  of  manipulation  the  non-linear  wave  equation  may  be 
approximated  in  terms  of  velocity  potential 

4«£i*  * [(*>*♦ 


(2.10) 


where 


It  should  be  noted  that  if  the  fluid  is  lossless  and 
Cp=Cg  then  (2.10)  reduces  to  a previously  known  non-linear 
wave  equation  C 9} 
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iiMflial 


(2.10a) 


In  order  to  express  the  approximate  non-linear  wave 
equation  in  terms  of  acoustic  pressure  and  particle  velocity, 
one  can  rearrange  the  Eq.(2.9)  in  terms  of  p and  u and 
combine  that  equation  with  (2. 10). The  result  is  a quadratically 
non-linear  wave  equation  £ 2] 


«<&)■ -7#  *<tn 


(2.11) 


If  it  chances  to  be  that  and  nearly  satisfy 

the  wave  equation,  C,  D ( )s0  ,then  on  the  RHS  of  Eq. 

(2.11)  * -^£1  and  (2.11)  becomes 

Further. if  it  happens  that  33  f 'P  f a 1 . 

3t*  ^ f.C}'  3t»  '•C.  / 

as  is  true  for  solutions  to  the  wave  equation  separated  in 
cartesian  coordinates,  then  [2] 


(2.12) 


(2.13) 


(Note  that  this  is  true  only  for  cartesian  coordinates.) 

As  it  is  stated  in  £3} "The  LHS  of  Eq.(2.13)  is  the 
classical , linear  wave  equation  pertinent  to  the  system 
under  study. The  RHS  can  be  interpreted  as  a forcing  function 
consisting  of  a three-dimensional  spatial  distribution  of 
phase-cohorent  sources. In  a second-order  perturbation  theory, 


H 


"this  forcing  function  is  obtained  from  the  classical  (first- 
order)  solution  of  the  acoustic  problem. The  second-order 
perturbation  solution  describes  the  nonlinearities  resulting 
from  the  self  interaction  of  the  classical  solution. Higher- 
order  perturbation  solutions  consider  the  interaction  of  the 
non-linear  solution  with  itself, and  the  forcing  function  is 
composed  of  products  of  both  classical  and  nonlinearly 
generated  terms . 

Thus, if  a system  is  driven  at  frequency  w,the  nonlinear 
term  in. .. "equation  (2.13)  "...  will  force  the  existence  of 
all  integer  multiples  nw  of  the  driving  frequency  and  the 
full  solution  must  contain  all  harmonics  of  the  input 
frequency. In  a closed  cavity, each  of  those  nonlinearly 
generated  waves  whose  frequency  lies  close  to  the  resonance 
frequency  of  a standing  wave  of  the  cavity  and  whose 
associated  spatial  function  matches  that  of  the  standing 
wave  can  be  strongly  excited. Just  how  strongly  will  depend 
on  the  quality  factor  Q for  the  particular  resonance  and 
the  difference  between  the  resonance  frequency  of  the 
standing  wave  and  the  harmonic  nw 

"Consider  two  limiting  cases. 

"(1)  If  the  forcing  function  does  not  have  its  frequency 
nw  close  to  wn,this  standing  wave  is  being  forced  at  a 
frequency  far  removed  from  its  resonance . Since  this  yields 
the  inequality 

I co Of  I » 

losses  can  be  ignored  in. . . "Eq. (2 .13) . 


.AJ* 
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'*(2)  if  nw*vwn,then  the  standing  wave  is  being  forced 

near  resonance , and  losses  must  be  retained  in. . . "Eq. (2 .13) . 

"The  value  of  Cn  can  be  determined  from  the  apparent  dimensions 

of  the  cavity  and  the  measured  resonance  frequency  wn. 

"The  losses  are  described  by  the  measured  Q of  the 

n 

resonance .This  means  that  the  linear-wave  equation  operator 

for  the  system  can  be  written  as 

2 __2  ^ a • r2  n2  n\v.  _JL 

CoD  +Tt£'  " C"D  ' Q»  9t  (2.15) 

"Comparison:'  of  cases  Cl)  and  (2)  reveals  that  the 
responce  of  the  cavity  when  nw^bw^  is  order  of  1/Q  compared 
to  that  when  nw/v/w^.Thus  for  the  high-Q  resonances  usually 
encountered  in  cavities  with  rigid  walls, the  components  in 
the  forcing  function  which  excite  standing  waves  far  from 
resonance  can  be  ignored  compared  to  those  components 
exciting  standing  waves  near  resonance.... 

"The  non-linear , coupled , transcendental  equation  appliciable 
to  this  problem  can  be  expressed  as 

Rnfsin|^n-8n>=N0Mof"5n  =°s9n  [1/2  RjVj{Iinj 

-£V)  Ri[sinj  (2.16) 


for  all  n>l.The  values  of  and  w must  be  determined  from 

^n  n 

the  infinitesimal-amplitude  behavior  of  the  cavity. The  Mach 
number  MQ  and  driving  frequency  w are  known  and  Nq  has  the  value 
Nq=1/2  for  a one-dimensional  standing  wave 
1/4  for  a two-dimensional  standing  wave 
1/8  for  a three-dimensional  standing  wave." 
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Rn  is  the  Fourier  coefficient  of  nth  harmonic  component, 
normalized  such  that  R^sl.  /n  is  the  phase  angle  of  the 
nth  harmonic  component , where  /^»0.,and  the  phase  angle  0n 
is  given  byf3] 


tan0  =-F_ 
n n 


(2.17) 


where  Fn=Qn[(nw ) 2-wn2]/ (nw) 2 =2Qn(SHz”n)  for  «1  (2.18) 

n wn 

Equation  (2.16)  can  be  solved  by  a method  of  successive 
approximations  on  a digital  computer . This  has  been  done  by 
[3]  and  [Si  and  both  decided  that  the  theoretical  model  can 
be  used  to  identify  the  modes  of  a non-ideal , rigid-walled 
cavity  provided  quantities  en  to  be  defined  later  are 
sufficiently  small. The  theoretical  model  in  its  present 
form  fails  to  account  for  the  excitation  of  modes  other 
than  those  belonging  to  the  family  of  the  driven  mode . 

This  excitation  was  observed  to  occur  only  in  the  case  of 
nearly  degenerate  modes.lt  is  believed  to  be  caused  by 
some  linear  coupling  mechanism  within  the  cavity. 

The  purpose  of  this  research  is  to  see  if  the  presence 
of  wall  irregularities  can  explain  how  non-family  members 
may  be  strongly  excited, and  to  present  an  example  to  support 
this  theory. 


3. DEFINITIONS  AND  NOTATIONS 


A.  FREQUENCY  PARAMETER 

The  frequency  parameter  is  a quantity  which  indicates 
the  position  of  the  driving  frequency  relative  to  the 
resonance  frequency , f^ ,of  driven  mode  in  terms  of  the 
of  the  driven  mode. The  frequency  parameter  is  defined  by 

Fl=2Ql(f"fl)/fl  (3.1) 

B.  STRENGTH  PARAMETER 

The  investigation  of  the  pressure  waveform  in  the  cavity 
required  the  calculation  of  the  strength  parameter  from  the 
observable  quantities . The  strength  parameter  is  defined  as 

STRPM=MQfQ1  (3.2) 

where 

VV^o'V’  <3'3) 

and  P^  is  the  peak  amplitude  of  p15the  pressure  distribution 
of  the  driven  mode. 

In  terms  of  observable  or  calculable  quantities  it  is 
reformulated  asT5j 

STRPM=>/7  V|Q1/(Sm/>C0)2  (3.4) 

where  V and  S^  are  the  rms  voltage  reading  and  microphone 
sensitivity  respectively  of  the  receiver  used  to  sense  the 
standing  wave. 


C.  e is  defined  to  indicate  the  position  of  w relative 
n n 


to  the  classical  harmonic  frequencies ,nw^ . 
w -nwi 

e = -2 i 

n nw^ 

and  one  can  relate  e with  F from  (2.18)  such  as 

n n 

Fn-2Vn 


(3.5) 


(3.6) 


D.  A pictorial  representation  of  which  will  be  useful 
throughout  the  development  is  given  in  Fig. 1. From  now  on 
three  subscripts  will  be  used  for  convenience,  i.e. 

Fn  becomes  Fnml.  . 


Q«ni 


I — ( ^nwt 
nw  / 


FIGURE  1 


F =fl_( 
nml  t nw 


wnml  . 2] n _ , . 

nw  -»^nml"  ° ^nml 


wnml  Sin«nml'(1+F2  )l/2‘bnml 

nml 


(3.7) 


(3.8) 
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4.  THEORETICAL  DEVELOPMENT 


A.  CAVITY  DESCRIPTION 

Assume  a perfectly  rigid-walled  rectangular  cavity  which 
has  one  wall  perturbed  such  that  the  cavity  dimensions  are 
Lx£l+ € f (y ,z)J  ,Ly  and  Lz  as  shown  in  Fig. 2 below. Also  assume 
the  perturbation  on  the  boundary  is  very  small  compared  to 
the  cavity  dimensions,  |6f(y,z)l«  1. 


FIGURE  2 

The  cavity  is  to  be  excited  by  a source  near  the  origin 
in  such  a way  that  the  (N,M,L)  mode  is  driven  at  a frequency 


close  to  its  resonance  frequency. 


B.  THE  PERTURBED  BOUNDARY 

For  a rigid-walled  rectangular  cavity  with  ideal  boundaries 
(€=0),the  pressure  pn  obtained  from  the  linear  wave  equation 


with  losses 

□l2  Po*° 

is  subject  to  the  following  conditions, 
Vpn  .n=0  at  x=Q , L 


(4.1) 


Vpg  ."n=0  at  x=0,Lx  (4.2) 

y=o,Ly 

2=o, lz 

where  "n  is  the  local  normal  to  the  ideal  boundary. The  solution 
for  pn  in  terms  of  Mach  number  is  given  by [4] 


„=Mn  cosk  x cosk  y cosk  z cos(wt+0) 

/>0c02  0 x y z (4.3) 

=Mq  (N,M,L/w,0) 

where  k’s  are  given  in  Eq.(2.2),  and 

w=c  (k  2+k  2+k  2)1/2  (4.4] 

p x y z 

If  the  cavity  has  perturbed  walls, the  solution  will  be 
in  terms  of  a summation  of  the  classical  linear  solution 
for  ideal  boundaries  plus  perturbation  correction  terms 


(4.4) 


due  to  the  irregular  boundary: 

p=Pn+Cp'+«2p' ’+ 


(4.5) 


Since  the  magnitude  of  the  boundary  perturbation  is  kept 
small, second  and  higher  terms  in  6.  can  be  considered 


insignificant , so  that 

p=p0+€p'  (to  first  order) 
and  p must  satisfy  the  following  conditions, 

□ 2p  = 0 


(4.6) 


(4.7) 
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(4.8) 


Vp»^=0  at  x=0  ,LxCl+^^(y  >z\l  (4.8) 

y=o,Ly 

z=0, Lz 

where  n,the  local  normal  to  the  real  surf ace, is  obtained  by 
taking  the  gradient  of  the  equation  for  the  boundary .given 
by  [2] 

"n=  V^x-Lx[l+€f  (y,z)jj  (4.10) 

Thus, to  the  first  order  in  6; 

t=i- CL  J-cl  34^  z (4.11) 

x ^ J x fa 

and  when  Eq.(4.11)  is  used  in  (4.8)  the  result  is 

r 3f>  ifai)  .ti, yl  -o  <4-12) 
L^T  x 7)2  ^2  ^ 2 ’oi  Jx  = LxC»4-£f(^.I)] 

A Taylor  series  expansion [4]  for  p evaluated  at  the  real 
boundary  L [l+6f(y,z)]  produces 

-ol  * f I +73r^|  €.Lxp(y.z> 

^‘x*LxCi+€f(s.»)3  **Lx 

+ C€Lxf  <**)]  + (4’13) 

2 ***  'XaLx 


(4.11) 


(4.12) 


(4.13) 


Substituting  Eq.(4.6)  into  RHS  of  ( 4 . 13 ) , taking  the  partial 
derivative  with  respect  to  x on  both  sides  and  keeping  the 
first-order  terms  in  6.  , yields  , x 

2Z\  + I -♦* cLxfc>,i;+"' 

*3*  lx«UD+Cf<9'*>]  ^ '**1*  Xcl*  *«L*  (4.14) 


Taking  the  partial  derivatives  with  respect  to  y and  z and 
using  exactly  the  same  procedure  gives 

2L  .2*1  + e|£l|  tL.fW.is) 

13 '**L>  3 ,'L* 
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In  order  to  find  the  contribution  to  p'  from  each  of  the 
terms, it  is  stated[2}that  for  a cavity  forced  by  a dynamic 

boundary  condition  at  a boundary 

2 2 
— * »-r 


c-’Dl"p,s0 


(4.20) 


and  the  dynamic  boundary  condition  from  the  (m,l)th  term  is 

c(wt+0) 


Dx 

then 


cosier  e 


L*  Lj  LX 


(4.21) 

t(wt+0+(fSmO 

e 

(4.22) 

(4.23) 

and  Snml  is  given  by  Eq.(3.8). 

Applying  this  solution  to  Eq.(4.18)  gives  the  complete 
solution  for  the  first-order  perturbation 

OQ 

9 9 

f * Z._Tml 

® - (4.24) 


m=o 

LsO 


/~WTiT  ^ (n.m.l  / va/»  © -t  ^nrnL  ) 

(^.rL*  n»o 
' Co  / 

. wzO 

and  combining  (4.24)  with  (4.6)  yields  the  total  acoustic 
pressure  in  the  cavity. 


<22- 


f = fN,M,L/W,8)--lr-  (n.m.l 

fe)L‘r.  (4.25) 

1*0 


If  this  is  near  a resonance.w^w  , ,then  this  term  will 

nml 

dominate  the  summation  and  all  other  non-degenerate  terms 
can  be  omitted . Consequently , the  Eq.(4.25)  becomes 

? = e-Mw)  (4  26) 

Cc:)  l* 
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5.  SPECIFIC  EXAMPLES 


A.  CAVITY  WITH  STEPPED  PERTURBATION 

Assume  that  the  rigid-walled  rectangular  cavity  given  in  Fig. 2 
is  perturbed  as  shewn  in  Fig. 3 belcw,and  also  assume  that  the  cavity 
is  driven  in  the  (0,1,0)  mode  resulting  (O,l,O/w,0)  standing  wave 
and  that  the  (0,2,0)  and  (1,0,0)  modes  are  (nearly)  degenerate. 


From  Fig. 3 the  equation  for  the  boundary  at  can  be  found. 

By  means  of  Eq.(4.8)  6 and  f(y,z)  can  be  written  as 

e » Vl* 

81  L'<*<L 

Since  the  (0,2,0)  and  (1,0,0)  modes  are  degenerate  the  emphasis 


(5A.1) 

(5A.2) 
(5A. 3) 


of  this  development  will  be  on  these  particular  modes. The  pressure 
distribution  of  (0,2,0)  node  is 

?oio  = Pi c 03  45  * cos C2vi * 9l ) 

or  (5A.4) 

-Poio  - ?!  (0,2,0/2W,62) 

where  is  the  amplitude  of  (0,2,0)  mode. 
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Utilizing  the  theory  developed  in  the  proceeding  sections  and 
using  the  equations  (4.17)  through  (4.26),  the  first-order  perturbation 


correction  can  be  found 


w 

xsL, 


, r ctfoio  1 

tj-Jx.l,  (5A*5) 

_ Lx  | ^cin  2TT  ^ COS(2Wt  + 0*)J[£()|-r)“£(ll"  L)jj  ^ 


Eq.(5A.  5)  can  be  written  as  a Fourier  series 


- 2TTftL*  V am  COS-p)/  C0sC2W-t +02)  ( 5 A . 5 ) 
l*1*.  U 


Inversion  of  the  Eq. (5A. 5 ) and  (5A.6)  yields  the  Fourier  coefficients 


3m  ,i[afe(2lL‘)co5(^L')-si«(iaL)cos(^L)]  (! 

and 

3o  ,-L[*n(&L‘) -**(&>-)]  (5A  8) 

Recalling  Eq.(4,23)  and  (4. 24), first  order  perturbation  correction 


(5A.7) 


is  found  as  . * f9w±  + e,+  ff^«) 

Va-iiL&ii  (2)(-l)a,005,n<,r00  COS^ii  e (5A'9) 

L>  (g) L*  * 

nd 

/ rrp  4 —v  <-(2W”fc  + 02  -tUioo) 

€f  = i£^6a0-ijriQ.IOo  °Too  e (5A.10) 

* (cr) 


where 


Hence, the  total  pressure, associated  with  the  angular  frequencv  2w, 
in  the  cavity  is 

f * fol 0 + &?' 

- F2(o,2,O/2M/;02)+  ^3L6aoQloo5ir»0^  (i,0,o/2W 

(5A.  ID 


The  total  pressure  at  the  ndcrophone  position ,x=0  and  y=Ly  ,is 

„ C t(2^we»)  r , - i (lwk +9*  + *wo)  7 

-pj  ^ + £ "^C^oQioo  ^,rt<^oo3  ® J (5A.12) 

Define  B=f  3 and  after  a little  manipulation  and  use  of  trigonometric 
identities, (5A.  12)  becomes 


'PLi*postt«r  P2|0+Bcoscroo)cos(zwUea).(6s,v,ar00)  sm(2wt+e,)j 

and  the  amplitude  of  the  total  pressure  in  the  cavity  is 

P Lic.p«it«n  ’ P2x/fl+fic««ro.)l+(B5mOr«,)1 


(5A.13) 


C5A.14) 


Eq.(5A.14)  is  the  corrected  value  of  the  amplitude  of  the  second 
harmonic  of  the  driving  mode , obtained  bv  Ea. (2. 16 ), because  of  the 
boundary  irregularity  given  in  Fig. 3. 

Now, it  is  desired  to  express  sin<J£#  in  terms  of  the  frequency 
parameter  of  the  driving  mode, (0,1,0). With  the  help  of  Fig. l.sinO^, 
can  be  written  as 

1 t — VU  ( 5 A . 15 ) 


If  f approaches  to  zero  then  CJ^Q  approaches  to  IT  ,and  if  f 
approaches  to  infinity  then  (7^0  is  close  to  zero. In  these  same 
limits  cosGTT  goes  to  -1  and  +1  respectively.  Hence 


For  2f  f 

sinC^g  = 


,(5A.15)  becomes 

i 


C5A.17) 


U+  UQioo  itf!22  Y V/a” 

L rio«  J 

Recalling  Eq.(3.1)  and  (3.5),F100  and  e1Qg  can  be  written  in  the  form  of 

(5A.18) 


and 


fiio  = *GU.o 

foio 

efioo  - 2ft)»o 
too  = — 1 


2-Poia 

Eq.(5A.19)  can  be  solved  for  f. 


100 


(5A.19) 

and  this  substituted  into  (5A.17) 


2Qm  = 2Qo..  5™. 

TIOO  QOIO  foio  C ' ■+  filOO  ) 


Qoio  foioC'  + fiioo) 

Use  of  l/(l+e10g)=l-e10g  and  little  manipulation  reveals 


(5A.20) 


S'*  Gii so  = 


^‘Ool  (l- e»oo)J 


(5A.21) 


As  a result , equations  (5A.14) ,(5A.16) ,(5A.21)  are  the  final  amplitude 
correction  of  the  second  harmonic  thi  driving  mode  obtained  bv  Ea . ( 2 . 16 ) . 
A computer  program  for  this  was  developed  bv  author  and  is  given  in 
appendix  A. 


B.  LINEARLY  PERTURBED  CAVITY 


Using  the  sane  assumptions  in  section  A, assume  that  the  rigid-walled 
rectangular  cavity  is  perturbed  linearly  as  shewn  in  Fig. 4 be lew. 


The  equation  for  this  perturbation  is 


X = L„  [■+•£•  (1--^  y )] 

(5B.1) 

Hence, 

(5B.2) 

and 

( 5 B . 3 ) 

Applying  the  same  procedure  as  in  section  A, the  ^irst-order  perturbation 
correction  and  the  total  acoustic  pressure  associated  with  angular  frequency 
2w  in  that  particular  cavity  can  be  found 

-Si  I * L-*  ? [- sin cos  (2Wfc +9*)]  (--£-)] 

(5B.4) 

= +i***st  C0s(2vt/i+9z) 

Li  ^ 


Eq. (5B.4)  can  be  written  as  a Fourier  series  and  the  Fourier  coefficient, 
am,is  found  bv  an  integration  procedure  evaluated  in  the  interval 
0 to  Ly  .The  result  is 


j_  f cos(2-*QTT  Cosq4m)TT  l 

TT  l (2 -ns ) (2  + m)  (2-iw) 


C24m)  J 


, 


So  s 0.0 

3l  a 0 0 
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(5B.5) 


Recalling  the  Eq,  (_4 , 24  l,the  first-order  perturbation 


correction  for  CQ.,2,Q1  mode  is 

f (2)C-i)lSm,(i.o,o/2vi,e1*o-um)  (5B.6) 

v (if;  l« 

and  the  total  acoustic  pressure  associated  with  angular 
frequency  2w  in  the  cavity  becomes 

since  aQ=0.0 

? = (0  .2,0/2<a>,9z)  ( 5B . 8 ) 

According  to  the  calculation  developed  above  there  is 
no  need  to  make  a first-order  perturbation  correction  to  the 
C0,2,Q1  mode  in  the  cavity  shown  in  Fig. 4,  As  a result, the 
pressure  distribution  is  equal  to  the  second  harmonic  of  the 
driven  mode, since  aQ=Q.Q  and  this  yields  -p'sQ.Q 

I 


f s (o»2,o/2uJ,0j ) 


'Sioo  (liO*O/2W,0j+OJoo)  ( 5 B . 7 ) 


C.  CAVITY  WITH  WEDGED  PERTURBATION 


Under  the  same  assumptions  made  in  sections  A and  B, 
assume  that  the  cavity  is  perturbed  as  shown  in  Fig. 5 below 


Applying  exactly  the  same  procedure  followed  in  section  A 
the  total  pressure  amplitude  in  the  cavity (in  first-order 
perturbation)  is 


P2VO+  Bcoso£0)  + (Bsin^oo) 


sin0J£o  and  costfj^  are  given  by  Eq.(5A.21)  and  C5A.16)  respectively. 

The  theoretical  predictions  of  these  specific  examples 
were  examined  with  series  of  experiments  developed  by  £8]  . 

The  further  discussions  about  these  will  be  given  in  the  next 


6.  RESULTS 


In  this  section  the  theoretical  predictions  performed 
in  sections  5A,  5B  and  5C  will  be  compared  to  the  experimental 
results  obtained  by  C83. 

The  information  on  the  empirical  losses  and  resonance 
frequencies  is  contained  in  the  Q's  and  e's. These  are  the 
values  used  in  the  computer  program  to  predict  the  harmonic 
distortion  on  the  basis  of  Eq.(2.16)  and  is  plotted  as  thin 
solid  curves. The  results  of  including  the  first-order 
perturbation  correction  are  plotted  as  thick  solid  curves 
for  each  specific  example. The  theoretical  curves  in  figures 
6 through  16  were  plotted  along  with  the  experimentally- 
measured  values  for  the  cavity  configurations  shown  on  top 
of  each  figure. The  theoretically-predicted  values  were 
generated  for  frequency-parameter  intervals  of  0.2,  and 
the  experimentally-measured  values  were  plotted  as  square- 
blocks.  Data  were  taken,  and  theoretical  predictions  made, 
for  different  strength  parameters  for  the  (0,n,0)  mode 
associated  with  different  cavity  configurations.lt  is 
important  to  note  that  the  n=2  distortion  peaks  when  the 
system  is  driven  at  this  frequency .That  is,  when  the  driving 
frequency  w is  equal  to  (1/2)  w2 , there  is  maximum  content 
of  Pj.The  point  where  this  occurs  for  each  P /P^  curve  is 
indicated  by^arrow  with  the  label  of  F^g.At  that  point 
the  value  of  frequency  parameter  is 

F020  = 2 -010  e020 

(The  same  thing  could  be  done,  of  course,  for  any  member 
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of  the  (0,n,0)  family).  The  arrow  labelled  as  F100  indicates 
the  position  where  the  nearly  degenarate  (1,0,0)  mode  is 
resonant,  and  the  value  of  F10q  is 

F100=2  Qoio  eioo 

The  theoretical  predictions  made  in  section  5B  were 
compared  to  the  experimental  results  as  seen  in  Fig. 7. 

When  Fig. 7 is  compared  to  Fig. 6,  the  unperturbed  cavity, 
it  is  clearly  seen  that  the  theory  and  experiment  are 
excellently  in  agreement. 

For  a wedged  perturbation,  the  theory  predicts  the 
frequency  of  the  second  harmonic  at  which  the  effect  of  the 
perturbation  occurs  as  seen  in  Fig. 8. The  predicted  magnitude 
of  the  perturbation  effect  for  this  configuration  is  in  good 
agreement  with  the  experiment . The  anomolous  behavior  of  the 
third  harmonic  in  Fig. 8 is  unexplained. 

For  stepped  perturbation,  when  the  cavity  is  perturbed, 
but  the  geometry  leads  to  no  predicted  correction  as  seen  in 
Fig. 9 or  leads  to  predicted  correction  less  than  about  0.02 
as  in  Fig. 12  or  less  than  about  0.05  as  in  Fig. 14,  then  it 
was  observed  that  there  was  very  little  or  no  effect  from 
the  (1,0,0)  mode .Agreement  for  these  cases  is  good  except 
for  the  region  lying  between  frequency  parameter  4 and  9 
in  Fig. 9. What  happened  in  that  region  is  also  unexplained, 
but  it  was  observed  one  time  only. When  the  amount  of 
perturbation  correction  is  increased  the  theory  predicts 
effects  larger  than  experimentally  observed . However , the 
effect  of  the  perturbation  appears  at  the  right  frequency 
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r- 


! 


parameter  as  is  seen  in  Fig.’s  10,11,13,15  and  16. 

Choosing  the  shim  position,  length  and  magnitude  is  very 
important  as  well  as  is  choice  of  the  strength  parameter. 

For  the  shims,  A =0.04  and  0.25  inches  for  stepped  and  wedged 
perturbations  respectively .The  effect  of  strength  parameter 
can  be  seen  in  Fig.'s  15  and  16. The  experimental  data  associated 
with  the  third  harmonic  in  Fig. 15  were  believed  to  come  from 
harmonic  distortion  in  the  piston  motion. 


’ 
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7.  CONCLUSION 


. JiilLUULWT" 


Non-linear  theory  has  been  applied  to  standing  waves  in  a 
rigid-walled  rectangular  cavity  with  a perturbed  boundary 
in  order  to  find  one  possible  mechanism  for  the  excitation 
of  a standing  wave  other  than  those  belonging  the  family  of  the 
driven  mode . It  was  observed  that  such  an  excitation  exists 
if  the  boundary  perturbation  and  the  dimensions  of  the  cavity 
are  favorably  chosen. 

It  appears  that  the  present  theoretical  model  succesfully 
predicts  the  major  features  of  harmonic  content  for  finite- 
amplitude  standing  waves  in  the  cavity  when  the  geometry  leads 
to  no  perturbation  correction  (Fig. 7 and  9). When  the  magnitude 
of  the  perturbation  is  increased  the  predicted  features  were 
larger  than  experimentally  observed . Second-order  perturbation 
corrections  may  be  needed  to  account  for  these  discrepancies. 

. 

^ 1 
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APPENDIX  A 


The  original  computer  program  for  Eq.(2.16)  was  prepared 
by  Coppens  in  1973,  and  author  made  some  extentions  to  that 
program  so  that  it  would  (1)  calculate  the  perturbation  cor- 
rection and  (2)  present  the  results  graphically.  The  program 
calculates  the  relative  amplitudes  and  phase  angles  of  stan- 
ding waves  in  cavities  keeping  the  strength  parameter  constant 
and  changing  the  frequency  parameter  to  generate  response 
curves  showing  the  amplitudes  of  the  nonlinearly  excited  stan- 
ding wav*s  as  function  of  the  frequency  parameter  . It  also 
calculates  the  perturbation  correction  according  to  Eq.(4.24) 
and  then  finds  the  total  relative  pressure  amplitude  using 
Eq.(4.26).  The  program  also  draws  the  graph  of  the  relative 
pressure  amplitudes  of  the  ideal  cavity, total  relative  pressure 
amplitude  of  the  perturbed  cavity  and  the  experimental  data 
on  a 3 cycle  semilog  paper. The  Versatec  Graphics  Plotting 
Manual  ClQ]  was  used  for  the  graphical  processes  on  the  IBM 
360  of  the  Randolph  Church  Computer  Center, Naval  Postgraduate 
School . 


SOME  USEFUL  INFORMATION  ABOUT  COMPUTER  PROGRAM 
* Quantities  marked  with  (*)  must  be  controlled  or  changed 
for  each  run. 

*K0N  The  number  of  iterations  throughout  the  region 

of  interest. For  this  program  the  iterations  are 
performed  with  0.2  intervals. 

*NCUR  The  number  of  the  experimental  curves  to  be 

drawn  + 1 

*NDAT  The  number  of  experimental  data  in  the  region 

of  interest 

BUR(I,J)  The  array  that  stores  the  experimental  data 
*XL  The  length  of  the  cavity  in  the  x-direction 

*YL  The  length  of  the  cavity  in  the  y-direction 

*DELTA  The  magnitude  of  the  perturbation 

*STRPM  Strength  parameter 

*FREQ  Frequency  parameter  stored  in  ATA(I,1)  and 

ZERCI ,1) . Input  as  the  maximum  value  of  FREQ 
in  the  region  of  interest 

XDAT(JET)  The  x coordinate 

YDAT(JET)  Value  of  the  curve  f(x)  for  XDAT 

ATA( I ,N) ,N^lThe  array  that  stores  the  logarithmic  value 
of  the  pressure  amplitudes  of  the  harmonics 
of  the  driving  mode 

ZER( I ,N) ,N^lThe  array  that  stores  the  linear  value  of  the 
pressure  amplitudes  of  the  harmonics  of  the 
driving  mode 

*Q(I)  Quality  factors  of  driving  mode  and  harmonics 

of  it 

*E(I)  e's  value  of  driving  mode  and  harmonics  of  it 
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*Q100 

*E100 

*XDAT( JET+1) 

*YDAT( JET+1) 

XDAT( JET+2 ) 

YDAT( JET+2 ) . 
HUM 

*B 


1 ’ 


Quality  factor  of  (1,0,0)  mode 
e-value  of  (1,0,0)  mode 

Integer  value  of  left-hand  corner  on  the  x-axis. 

It  must  have  the  same  value  as  the  7th  argument 
of  subroutine  CALL  AXIS  for  x-axis. 

Integer  value  of  left-hand  corner  on  the  y-axis . 

It  must  have  the  same  value  as  the  7th  argument 
of  subroutine  CALL  AXIS  for  y-axis. 

Increment  value  of  x and  y for  scalling  purposes 

The  linear  value  of  the  total  acoustic  pressure 
amplitude  associated  with  angular  frequency  2w. 

JLqi00  sin(J^00(l/4)  L aQ  for  stepped  perturbation 

Q100  sin(j^"00(a0/2 ) for  wedged  perturbation 

and  aQ  is  the  0th  Fourier  coefficient. 

I 1 1 
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